Decoherence of cold atomic gases in magnetic micro-traps 
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We derive a model to describe decoherence of atomic clouds in atom-chip traps taking the excited 
states of the trapping potential into account. We use this model to investigate decoherence for a 
single trapping well and for a pair of trapping wells that form the two arms of an atom interferometer. 
Including the discrete spectrum of the trapping potential gives rise to a decoherence mechanism with 
a decoherence rate F that scales like F ~ 1/ro with the distance ro from the trap minimum to the 
wire. 



I. INTRODUCTION 

Cold atomic gases form an ideal system to test fun- 
damental quantum mechanical predictions. Progress in 
laser cooling made it possible to achieve previously in- 
accessible low temperatures in the nK-range (see, e.g., 
the Nobel lectures 1998 (jy). One of the most exciting 
consequences of this development has been the creation 
of Bose-Einstein condensates (BECs) 0, Q . These have 
been manipulated by means of laser traps in various man- 
ners, e.g. vortices have been created or collisions of two 
BECs have been studied 0, 0, An interesting link 
to solid-state phenomena has been established by creat- 
ing optical lattices, in which a Mott transition has been 
theoretically predicted and observed Q. 

Recently, proposals to trap cold atomic gases using 
microfabricated structures have been realized exper- 
imentally [H in 113 on silicon substrates, so-called 
atom chips. These systems combine the quantum me- 
chanical testing ground of quantum gases with the great 
versatility in trapping geometries offered by the micro- 
fabrication process. Micro-fabricated traps made it pos- 
sible to split clouds of cold atomic gases in a beam split- 
ter geometry , to transport wavepackets along a con- 
veyer belt structu re [l4l and to accumulate atomic clouds 
in a storage ring [l^- Moreover, a BEC has been suc- 
cessfully transferred into a micro-trap and transported 
along a waveguide created by a current carrying micro- 
structure fabricated onto a chip 0, 0, [iS^ J^J . Finally 
several suggestions to integrate an atom interferometer 
for cold gases onto an atom chip have been put forward 
[20,21,22]. 

The high magnetic-field gradients in atom-chip traps 
provide a strongly confined motion of the atomic quan- 
tum gases along the micro-structured wires. The atomic 
cloud will be situated in close vicinity of the chip sur- 
face. As a consequence, there will be interactions be- 
tween the substrate and the trapped atomic cloud, and 
the cold gas can no longer be considered to be an iso- 
lated system. Recent experiments reported a fragmenta- 
tion of cold atomic clouds or BECs in a wire waveguide 
[T^ I23L |3j us on reducing the distance between the 
wavepackets and the chip surface, showing that atomic 
gases in wire traps are very sensitive to its environment. 
Experimental ^2&\ and theoretical works showed that 



there are losses of trapped atoms due to spin flips induced 
by the strong thermal gradient between the atom chip, 
held at room temperature, and the cold atomic cloud. 
Moreover, Refs. 28, 29] studied the influence of magnetic 
near-fields and current noise on an atomic wavepacket. 
Consequences for the spatial decoherence of the atomic 
wavepacket were discussed under the assumption that the 
transversal states of the one dimensional-waveguide are 
frozen out. 

For wires of small width and height, as used for atom- 
chip traps, the current fluctuations will be directed along 
the wire. Consequently, the magnetic field fiuctuations 
generated by these current fiuctuations, are perpendic- 
ular to the wire direction, inducing transitions between 
different transverse trapping states. 

In this paper we will discuss the influence of cur- 
rent fluctuations in micro-structured conductors used in 
atom-chip traps. We derive an equation for the atomic 
density matrix, that describes the decoherence and equi- 
libration effects in atomic clouds taking transitions be- 
tween different transversal trap states into account. We 
study decoherence of an atomic state in a single waveg- 
uide as well as in a system of two parallel waveguides. De- 
coherence effects in a pair of one-dimensional (Id) waveg- 
uides are of particularly high interest as this setup forms 
the basic building block for an on-chip atom interferom- 
eter iMmm. 

Our main results can be summarized as follows. Con- 
sidering current fluctuations along the wire, we show that 
spatial decoherence along the guiding axis of the Id- 
waveguide occurs only by processes including transitions 
between different transversal states. The decoherence 
rate obtained scales with the wire-to-trap distance rg as 
F l/''o- Applying our model to a double waveguide 
shows that correlations among the magnetic-field fiuctu- 
ations in the left and right arm of the double waveguide 
are of minor importance. The change of the decoherence 
rate under the variation of the distance between the wires 
is dominated by the geometric rearrangement of the trap 
minima. 

To arrive at these conclusions we proceed as follows. 
Section^will give a derivation of the kinetic equation de- 
scribing the time evolution of an atomic wavepacket in an 
array of N parallel Id waveguides subject to fiuctuations 
in the trapping potential. We will then use this equa- 
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(a) 




FIG. 1: Setup of the atom chip showing directions of the 
wires and the magnetic bias fields needed to form the trapping 
potential, (a) Single-wire trap, (b) Double-wire trap. 

tion of motion in Section UTTI to study the specific cases 
of a single Id waveguide and a pair of Id waveguides. 
Spatial decoherence and equilibration effects in the sin- 
gle and double waveguide will be discussed. Finally we 
summarise and give our conclusions in Section [IVI 
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FIG. 2: Contour-plot of the magnetic field of the single-wire 
trap for B^f^^ = lOG, B^H^ = 5G, and I = 0.1 A. The up- 
per (right) plot shows a horizontal (vertical) cut through the 
potential minimum. The dashed lines in the upper and the 
right plot show the harmonic approximation to the trapping 
field. 



II. EQUATION FOR THE DENSITY MATRIX 



Neumann equation 



The trapping potential of a micro-structured chip is 
produced by the superposition of a homogeneous mag- 
netic field and the magnetic field induced by the current 
in the conductors of the micro-structure. Atoms in a 
low- field seeking hyperfine state \S) will be trapped in 
the magnetic field minimum |3fl| | . The interaction of the 
atom with the magnetic field is 



V^(x)=-(5|/x|5)B(x), 



(1) 



assuming that the Larmor precession of the trapped atom 
is much faster than the trap frequency and, hence, the 
magnetic moment of the atom fi can be replaced by its 
mean value (Sj/xIS*). 

Many different wire-field configurations will lead to an 
atom trap . We will concentrate on systems in 

which the trapping field is produced by an array of par- 
allel wires, see Fig. ^for a single or double-wire trap. A 
homogeneous bias field is applied parallel to the surface 
on which the wires are mounted. An example for the re- 
sulting field distribution for the single- wire trap is shown 
in Fig. El 



A. Stochastic equation for the density matrix 

The quantum mechanical evolution of a cold atomic 
cloud is described by a density matrix. The time evolu- 
tion of the density matrix p(x, x',i) is given by the von 



(2) 



where H is the Hamiltonian of an atom in the trapping 
potential: 



H=-^ + Vt{ri_) + 5V{^,t) 
2m 



(3) 



Here, Vt{r±) is the confining potential which we will as- 
sume to be constant along the trap, and denotes the 
coordinate perpendicular to the direction of the current- 
carrying wire or the wave guide. The last term 5T^(x, t) is 
a random fluctuation term induced by the current noise. 

We will formulate the problem for a system of N par- 
allel quasi-ld magnetic traps generated by a set of M 
parallel wires on the chip (in general, N ^ M since some 
of the magnetic-field minima may merge). The number of 
parallel trapping wells is included in the structure of the 
confining potential Vt{v±). In all further calculations we 
will assume that the surface of the atom chip is in the x-z 
plane and that the atoms are trapped in the half space 
y > Q above the chip. The wires on the atom chip needed 
to form the trapping field, are assumed to be aligned in 
z-direction. 

In position representation for the density matrix the 
von Neumann equation reads 



^ d 



SV{:>i,t) - (5\/(x',i)jp(x,x',i) . 
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To derive a quasi-ld expression for Eq. we expand the 
density matrix in eigenmodes of the transverse potential 
Xn{r±): 

p(x,x',t) =^^Xni^l-)Xmir'j^)Pnmiz,z',t) . (5) 



Here, the channel index n labels the transverse states 
of the trapping potential. The transverse wavefunctions 
Xn(r±) are chosen mutually orthogonal and are eigen- 
functions of the transverse part of the Hamiltonian in 
Eq. (O in the sense that 



we will hereafter abbreviate this part by {ifidt — Hik)- 
The terms on the right-hand side of Eq. {Tj) contain the 
influence of the potential fluctuations. 

The fluctuations of the confining potential are de- 
scribed by the matrix elements 5Snk{z,t) which imply 
transitions between different discrete transverse energy 
levels induced by the fluctuating potential. The influence 
of the fluctuating potential onto the longitudinal motion 
of the atomic cloud is included in the z dependence of 

6Snk{z, t). 



Xn{rj_) = EnXnir±) ■ (6) 



The decomposition jSJ in transverse and longitudinal 
components of the density matrix is now inserted in the 
von Neumann equation, Eq. Making use of the or- 
thogonality and the completeness of the transverse states 
Xn(r_L) we obtain a one-dimensional equation for the evo- 
lution of the density matrix: 





















2m 





pik{z,z',t) (7) 

= ^ [SSlniz, t)pnkiz, z', t) - SSnkiz', t)pin{z, z', t)] . 



Here the abbreviation AEik = Ei — Ek was introduced 
for the difference of transverse energy levels, and the fluc- 
tuations dSnk{z) are defined as 

SSnk{z,t) = J dv^ Xn{ri_)5VMxl{vi_) . (8) 

The left-hand side of the reduced von Neumann equa- 
tion, Eq. |(2J), describes the evolution of an atom wave- 
function in the non-fluctuating trapping potential and 



B. Averaged equation of motion for the density 
matrix 



In this section we derive the equation of motion de- 
scribing the evolution of the density matrix {p) averaged 
over the potential fluctuations [s^ 01 ■ To keep the ex- 
pressions compact we rewrite Eq. {TJ as 



[ihdt -H^j p = SSp, (9) 
where the spatial coordinates x, x', time t and all indices 
have been suppressed. The components of the quantities 
H and SS are defined as 

Etikij = ElikSiiSkj , (10) 
SSikij = SSii{z,t)Skj ~ SSjkiz' ,t)Sii . (11) 

The product in ^ has the meaning Ap = AikijPij- 
The stochastic equation Q for the density matrix is av- 
eraged over the fluctuating potential using the standard 
cumulant expansion |34j . Writing the time arguments 
again we obtain 



[thdt-H) {p{t)) ^ {SS{t)){p{t)) -'-j\t'{{5S{t)e~i^''5S{t-t')))ei^''{p{t)) . (12) 

Here, the brackets ( • ) denote the averaging over all realisations of the potential fluctuations and the double brackets 
((•)) denote the second cumulant. We can take the mean value of the potential fluctuation {5V) to zero, since the 
static potential is already included in the Hamiltonian on the left hand side of Eq. (|12|l and, hence, the flrst term on 
the right-hand side of Eq. p2|l vanishes. 

Reinserting the explicit expressions for iJ, 5S leads to the desired equation of motion for the density matrix 

i oo oo 

(indt-Hik){pik{z,z',t)) =-^Y^ J dr J dS J di'K^j{z^ 2,z' ~ Z',t-T) (13) 

-oo -oo 

X [{SSli{z,t)SSi„i{z,T))6kj6jn + {SSjk{z\t)SSnj{z\T))SiiSi,n 

-{6Sli{z,t)6Snj{z ,T))5im5jk - {SSjkiz ,t)5Simiz,T))SjnSli] (prnn{z,Z,T)) . 

I 

The kernel Kij{z — z^z' — z' , t) is the Fourier transform which can be explicitly evaluated, but this has no advan- 
of tage for our further discussion. 



K,,{q, q\t)^ exp ( -t— l^q^ - q' y - -AE,jt 



(14) 
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Equation H13() is the main result of this section. It de- 
scribes the evolution of the density matrix in a quasi- Id 
waveguide under the influence of an external noise source. 
It is valid for an arbitrary form of the transverse confining 
potential, thus, in particular, for single and double-wire 
traps. The main input is the external noise correlator 
and its effect on the trapped atoms. The noise correlator 
depends on the concrete wire configuration. Below we 
will derive a simplified form of the equation of motion 
((T^ under the assumption that the time scale charac- 
terizing the fluctuations is much shorter than the time 
scales of the atomic motion. 

C. Noise correlation function 

We will now derive the noise correlator for our specific 
system and use Eq. H13|) to study the dynamics. As we 
are interested in the coherence of atoms in an atom-chip 
trap, we will consider the current noise in the wires as 
the decoherence source. Fluctuations of the magnetic 
bias fields B^^^^^ and B^^f^^, needed to form the trapping 
potential, will be neglected. 

Using the approximation of a Id wire, the fluctuating 
current density can be written as 

j(x,i)-/(z,<)J(x)%)z, (15) 

where z is the unit vector in z-dircction. The fluctuations 
of the current density are defined by (5j(x, t) = j(x, f) — 
(j(x)). It is sufficient to know {61{z, t)SI{z' , t')) to obtain 
the full current density correlation function ((5j(x)(5j(x')). 
Note, that the average currents are already included in 
the static potential 14 (rj.)- 

The restriction of j to the z direction is a reasonable as- 
sumption since we consider micro-structured wires with 
small cross-section A = l^hn i-c wires with widths 
and heights Ih much smaller than the trap-to-wire dis- 
tance tq. Transversal current fluctuations lead to surface 
charging and thus to an electrical field which points in 
opposite direction to the current fluctuation. This sur- 
face charging effect will suppress fluctuations which are 
slow compared to lorc — crA/eolw Here, a is the con- 
ductivity of the wire, eg the (vacuum) dielectric constant, 
and A « lA is the screening length in the metal. This 
leads to RC-frequencies of llirc ~ lO^'^Hz for wire widths 
of Iw = 10/J,m and typical values for the conductivity in 
a metal. 

The characteristic time scale for the atomic motion in 
the trap is given by the frequency of the trapping po- 
tential uj « lO'^Hz. Thus, considering atomic traps with 
uj ujRC the current fluctuations can be taken along z 
as a direct consequence of the quasi-one-dimensionality 
of the wire. 

The current fluctuations are spatially uncorrelated as 
they have their origin in electron scattering processes. 
Hence, the correlator {SI{z)SI{z')) has the form [sslls^ 

{SI{z, t)6I{z', t')) = 4fcBTcff (z)crA 5{z - z')5c{t -t') . 

(16) 



Here, fee is the Boltzmann constant. The effective noise 
temperature is given by [s^ 

Teff(z) = j dEf{E, z)[l - f(E, z)] , (17) 

where f{E, z) is the energy- and space-dependent non- 
equilibrium distribution function. A flnite voltage across 
the wire induces a change in the velocity distribution 
of the electrons and the electrons are thus no longer in 
thermal equilibrium. Nevertheless, the deviation from 
an equilibrium distribution is small at room temperature 
due to the large number of inelastic scattering processes. 
The effective temperature T^fi accounts for possible non- 
equilibrium effects such as shot noise. However, contribu- 
tions of non-equilibrium effects to the noise strongly de- 
pend on the length L of the wire compared to the charac- 
teristic inelastic scattering lengths. E.g. strong electron- 
phonon scattering leads to an energy exchange between 
the lattice and the electrons. The non-equilibrium dis- 
tribution is 'cooled' to an equilibrium distribution at 
the phonon temperature. Thus, non-equilibrium noise 
sources, such as shot noise, are strongly suppressed for 
wires much longer than the electron-phonon scattering 
length Icp and the noise in the wire is essentially given 
by the equilibrium Nyquist noise 37, 38]. As the wire 
lengths used in present experiments are much longer than 
Zop, we Toff « 300K in all our calculations. 
Finally, 

is a representation of the delta function. The correla- 
tion time Tc is given by the time scale of the electronic 
scattering processes. 

D. Simplified equation of motion 

The dominating source of the current noise in the wires 
is due to the scattering of electrons with phonons, elec- 
trons and impurities. These scattering events are corre- 
lated on a time scale much shorter than the characteristic 
time scales of the atomic system. This separation of time 
scales allows us to simplify the equation of motion l|13|l . 

As a consequence of Eq. (|16|l , the correlation function 
of the fluctuating potential will be of the form 

{5VM5V{-x',t'))^5,{t-t'){5V{^)5V{^)) . (19) 

Using Eq. ^ we obtain 

{5S,^{z,t)5Sn,{z' ,t')) = 6,{t - t'){6S,^{z)5SnAz')) 

(20) 

for the fluctuations of the projected potential. 

We replace all the correlation functions in the equa- 
tion of motion (|13|l by the expression (|20|l . The time 
integration can be performed using the fact that the av- 
eraged density matrix (pir)) varies slowly on the correla- 
tion time scale Tc. Thus, {p{t)) can be evaluated at time 
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t and taken out of the time integral. Finally, taking Tc Performing the remaining spatial integrations over z and 
to zero, the Fourier transform of the kernel 1)14(1 leads to z', Eq. H13|) reduces to 
a product of delta functions in the spatial coordinates. 



{ihdt-Hlkj {pik{z,z' ,t)) ^ --^^[{SSln{z)SSnr,i{z)){prnkiz,z',t)) 

inn 

+ {SSrnk{z')SSn,niz')) {pin{z, z' , t)) -2 ((55/„i (z)(5S'„fc (z')) {pmn{z, z' , <))] 

I 



(21) 



The trapping potential has N minima, i.e., the channel 
index n can be written as n = {a,nx,ny), a = 1, ...,N. 
We will now assume that the minima are well-separated 
in iy-direction such that we can neglect all matrix ele- 
ments SSij with different trap labels a. Under this as- 
sumption, the first and second term on the right-hand 
side of Eq. H21|) depend only on a single trap label. In 
the third term, the trap labels may be different for SSim 
as compared to SSnk- Since current noise in one par- 
ticular wire generates fluctuations of the magnetic field 
in all trapping wells, there will be a correlation between 
the fluctuations in different traps. It is hence this third 
term in Eq. (|21|l which describes the correlation between 
potential fluctuations in different traps. 

Using Eq. it is now possible to describe decoher- 
ence induced by current noise in trapping geometries of 
one, two or more parallel wires [3ll'32^. In the following 
sections we will discuss decoherence in two specific trap- 
ping configurations: the single-wire trap in Section flll Al 
and the double-wire trap in Section IIII Bl which is of 
particular interest for atom interferometry experiments 



III. 



TRAPPING GEOMETRIES 



We will now consider two specific trap configurations, 
the single-wire trap and the double-wire trap. The dy- 
namics of the noise-averaged density matrix will be dis- 
cussed using Eq. H21|) which was derived in the previous 
section. We assume that the wire generating the poten- 
tial fluctuations is one-dimensional. This assumption is 
reasonable for distances of the trap minimum to the 
wire much larger than the wire width ly^ and wire height 
Ih- Typical length scales are rp ~ lO/xm to 1mm and 
wire widths and heights 1^ ~ 10/im to 50/^m, Ih ~ Ipm 

M 

To keep the notation short, we will suppress the brack- 
ets denoting the averaging. Thus, from now on, the den- 
sity matrix pij denotes the average density matrix. 



A. The single-wire trap 

We are now looking at a specific magnetic trapping 
field having only a single one-dimensional trapping well 
as shown in Fig.|21 The magnetic field is generated by the 
superposition of the magnetic field due to the current / in 
the single wire along the z-axis and a homogeneous bias 

(x) 

field parallel to the chip surface and perpendicular 

to the wire (see Fig. 01. We additionally include a ho- 

mogeneous bias field S^ias parallel to the trapping well. 
This longitudinal bias field is experimentally needed to 
avoid spin flips at the trap center. Thus the magnetic 
field has the form: 



B(x) 



Mo J 



1 



2tt a;2 



-y 

X 





B 



B 



(x) 
bias 





(22) 



The minimum of the trapping potential is located above 
the wire — Q and has a wire-to-trap distance of tq = 

yo=Mo//(27ri?(^J. 

We have now all the ingredients needed to calculate the 
fluctuation-correlator {5Sim{z)5Snj{z')) . Using Eqs. l|SJl, 
(fTK|l . and (fTB|l we obtain 



{5Sira{z)5Sn]{z')) ^ A^mn] J [z - z') 

The transition matrix elements are given by 



(23) 



A,, 



k-QTcaa-A 



( pLogpP'B B 



(x) 



B 



(z) 



X j dY±_ Xi(i"±)(y-yo)Xm(r_L) 
X / rfrYxn(rl)(y'-yo)x:W) 



(24) 



The spatial dependence of the noise correlator is given 
by 



— OO 

(25) 

A more detailed derivation of the correlation function is 
given in Appendix ^ The integration in Eq. 1241) can 
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be done explicitly using harmonic oscillator states for 
Xn(r_L), which is a good approximation as long as B^^!^^ 

(x) 

is of the order of B^^J^^ and the wire-to-trap distance ro is 
much larger than the transverse width w of the trapped 
state. Both requirements are usually well satisfied in ex- 
periments. The dashed line in Fig.|2]shows the harmonic 
approximation to the trapping potential. Using the har- 
monic approximation we can characterize the steepness 
of the trapping potential by its trap frequencies. It turns 
out that the two frequencies coincide, 

(26) 

i.e., the trap potential can be approximated by an 
isotropic 2d harmonic oscillator (2d-H0). 
After integration we obtain 

where 

Styniy = ^my + l (5j^,„l^ + l + Sty,-my-l , (28) 

and the indices riy denote the energy levels in a; and 
y direction of the 2d-H0. The prefactor is 

\ bias / 

and w = ^hjimiij) is the oscillator length of the har- 
monic potential. 

Before moving on to derive the equation of motion for 
the averaged density matrix, let us discuss some conse- 
quences of expression (|27f) . Inspection of the coefficients 
Aimnj in Eq. H27|l shows that there is no direct infiu- 
ence onto the longitudinal motion of the atom, since the 
matrix element Anjj vanishes. This result is not un- 
expected, because we consider only current fluctuations 
along the wire, which give rise to fluctuations in the trap- 
ping field along the transverse directions of the trap po- 
tential. In fact, only transitions among energy levels 
of the y-component of the 2d-H0 give a non-vanishing 
contribution. An explanation can be obtained by exam- 
ining the change of the trap minimum position under 
variation of the current in the wire. Changing the cur- 
rent by SI leaves the trap minimum in ^-direction un- 
changed at xo — 0, but shifts the y-trap minimum by 
Syo = I-IoSI/{2ttBI^^IJ 39]. Even though there is no di- 
rect coupling to the motion along the wire there can still 
be spatial decoherence in z-direction as we will see in 
the final result for p obtained from Eq. (|21|) . However, 
this requires transitions to neighbouring transverse en- 
ergy levels. 

Substituting Eqs. and in Eq. leads to 
an equation of motion for p for the single- wire trap con- 
figuration. Instead of discussing the general equation of 
motion of p we will only consider the case where the two 




FIG. 3: Spatial dependence of the decay rate fdec(C-) = 
^(0) - A(C-) for the diagonal elements pii{C-,t), Eq. 
The wire-to-trap distance is ro = lOO^m and the trap fre- 
quency is oj « 27r X lOkHz. The dashed line is the approxima- 
tion Eq. 1391 for Inset: Decay rate A(0) as function 
of ro. The parameters chosen in the plots correspond to *^Rb 
trapped at TeS = 300K in a magnetic trap with a gold wire 
of conductivity (tau = 4.54 • lO^fi-^m"^ and a cross-section 
of yl = 2.5/im x 5/im. The bias fields chosen are iJ^iia ~ ^0^1 
and B^fi, = 2G. 



lowest energy levels of the transverse motion are taken 
into account. The restriction to the lowest energy levels 
corresponds to the situation of the atomic cloud being 
mostly in the ground state of the trap and having negli- 
gible population of higher energy levels. This situation is 
realistic, if the energy spacing of the discrete transverse 
states is large compared to the kinetic energy. Neverthe- 
less, the result obtained from the two-level model pro- 
vides a reasonable estimate for the decoherence of an 
atomic cloud, even if many transverse levels are pop- 
ulated. Equation (|27|l shows that only next-neighbour 
transitions are allowed. As Aimnj is a product of two 
next-neighbour transitions there can only be contribu- 
tions of the next two neighbouring energy levels. Higher 
energy levels do only contribute to the decoherence by 
successive transitions which are of higher order in Aimnj 
and hence are negligible. 

We rewrite Eq. H21|) for the subspace of the two lowest 
eigenstates in a matrix equation for the density vector 

P = (poo,Pii,Pio,Poi) • (30) 

The indices of the averaged density matrix pik denote the 
transverse state I = {lx,ly) but we are only writing the ly 
component of the label as p can only couple to states with 
the same x-state (i.e. only transitions between different 
j/-states of the 2d-H0 are allowed) . The x-label is = 
for all states under consideration. Hence Eq. H21|l can 
be written as a matrix equation for the density matrix 
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vector Eq. ^ 

ihdt — H{z, z') p{z, z' , t) = —iTiA{z — z')p{z, z', t) . 

(31) 

The matrix A is defined as 



MO) 





A{0) 
^(0) 

-A(c-) 




where A(C_) = J(C-) and C- 



(32) 
The equa- 
tions for the diagonal elements, i.e., poo, pn, and for the 
off-diagonal elements, i.e. poi, Pio decouple. However, 
the decoupling is a consequence of the restriction to the 
two lowest energy levels and is not found in the general 
case. Yet, the decoupling allows to find an explicit so- 
lution for the time evolution of the diagonal elements. 
The matrix equation proves to be diagonal for the linear 
combinations p^ = poo ± Pii and introducing the new 
coordinates C,- ^ z — z' , ~ ^(z + z') we obtain the 
general solution 



where the decay is described by 



(33) 



r±(C-,t) - A{0) T-J dt'AiC- - ^t') . (34) 



The function is fixed by the initial conditions, i.e. by 
the density matrix at time t = : 

oo 

i?±(C-)- J rfC+e-^'=^+p±(C-,C+,i = 0). (35) 

— oo 

Note that F in Eq. H34|) is a function of the spatial vari- 
able (- and the wavevector k, and is in general not linear 
in the time argument t. Adding and subtracting the con- 
tributions p^ finally leads to the following expressions for 
the diagonal elements of the density matrix 



Poo(C-,C+,<) = 



dk 
2^ 



[p+{C-,C+,t) + pj:{C-,C+,t)]. 

(36) 

The result for pn can be obtained from (|36|l by replacing 
the plus sign between the terms by a minus sign. To get a 
feeling for the spatial correlations we trace out the center 
of mass coordinate C+ : 



Poo{C-,t) = / dC+ Poo (C- , C+ , ^ ) 

^ — OO 

^ g-[A(0)-A(C-)]t 

1 



2 

+ Pii 



(C-,0)(l-e-2^«-)*) 




FIG. 4: Time evolution of |poo(C-, C+, *)! after t = 2/A{0). 
At t — the wavepacket is a Gaussian wavepacket of spatial 
extent Wz ~ 20ro. The wire-to-trap distance is ro — 5/im 
and all other parameters correspond to those used in Fig. |3| 
A{0) ~ 0.5s~^ is given by Eq. H88|l . The density matrix shows 
damped oscillations which become more pronounced as ("+ 
increases. 



Equation 1)3 7|) shows that we can distinguish two decay 
mechanisms. There is an overall decay of the spatial off- 
diagonal elements with a rate rdcc(C-) — [^(0) — ^(C-)]- 
This decay only affects the density matrix for C_ ^ 0, 
thus suppressing the spatial coherence. The spatial corre- 
lation of the potential fluctuations can be read off Fig. O 
showing fdoc(C-)- We find for the potential fluctuations 
a correlation length of the order of the trap-to-chip 
surface distance ~ ''o- The correlation length must 
not be confused with the coherence length describing the 
distance over which the transport of an atom along the 
trapping well is coherent. This coherence length is de- 
scribed by the decoherence time and the speed of the 
moving wavepacket. The correlation length character- 
izes the distance over which the potential fluctuations are 
correlated. 

The second mechanism describes the equilibration of 
excited and ground state, quantified by the diagonal ele- 
ments poo(0,i), which occurs at a rate fpop(O) = 2^(0). 
Here equilibration means that, due to this mechanism, 
the probability to be in the ground state, i.e. poo(0,i), 
tends to 1/2. Of course, at the same time the probability 
pii(0,t) to be in the excited state approaches 1/2 at the 
same rate. Note, that for the spatial off-diagonal matrix 
elements the equilibration rate depends on 

We will now take a closer look onto the quantity A{0) 
for realistic trap parameters. Using J(0) = ^ and 

Eq. |2ni), Eq. (EHl) leads to 



2pB9F 



B 



(^) 

bias 



(37) 



(38) 

The inset of Fig. O plots A(0) over ro for reasonable trap 
parameters. Equation (|38|l shows that the decay rate is 
scaling on the wire-to-trap distance as 1 /rjj giving a rapid 
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FIG. 5: Decay of the absolute value jpoo(C-iC+ — 30ro,t)| 
and Ipoo(C-) C+ = Oi ^)l (inset). The wave packet is a Gaussian 
wavepacket of spatial extent Wz = 20ro at t = 0. The wire-to- 
trap distance is ro = 5/im. All other parameters are the same 
as used in Fig. A(0) « 0.5s"^ is given by Eq. lEHJ- The 
dashed line (t = 2/A{Q)) corresponds to cuts of Fig. 0|along 
for ^-|- = 30ro (inset C,+ = 0). Spatial correlations of the 
potential fluctuations are restricted to a narrow band around 
C,- = 0. The width of this band is on the order of ro. The 
density matrix shows damped oscillations which become more 
pronounced with increasing (,+■ The increase of |poo(C-i C+ — 
30ro ,t) \ for increasing time t is a consequence of the spreading 
of the wavepacket. 



increase of decoherence effects once the atomic cloud is 
brought close to the wire. 

As specific example we want to study the time evolu- 
tion of the full density matrix in the ground state. Fig- 
ure 01 and Fig. El show the time evolution of the absolute 
value of the density matrix element |poo(C-j C-i-j for ^ 
Gaussian wave packet of spatial extent Wz = 20ro and 
a wire-to-trap distance of rg = 5^m. Initially, at time 
t = all other elements of the density matrix are zero. 
The time evolution is calculated for the reduced subspace 
using Eq. (|33|l and Eq. (|36|l . For the spatial correlation 
we use the approximation (dashed line in Fig. j^J 



influence of the fc-dependent damping we assume that 
r± = A(0) T [^(C-) + (i{C,-,t)k]. Choosing F linear in k 
leads to a modulation of the density matrix by a factor 
proportional to cosh(iQ;o/3 + ai), describing damped os- 
cillations. The ckQ/i s-re real functions of C± and t. The 
oscillations arise only for non-vanishing (3. Decay rates 
which are /c-independent do not show oscillations. 

The decay described by F^ in Eq. H34() is however not 
linear, but includes higher powers in k. Hence the simple 
linear model F^ = A{Q) =F [^(C-) + describes 
the oscillations in Fig. 01 and Fig. [SI only qualitatively. 
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FIG. 6: Contour-plot of the magnetic field in the double-wire 



geometry for B'^^.s = 10^, B^f^ = 5G, and I = 0.3A. The 
wires are located at a; = ±100/im and j/ = 0. Two poten- 
tial minima are located roughly above the current- carrying 
wires. The upper (right) plot shows a horizontal (vertical) 
cut through the right potential minimum. The dashed lines 
in the upper and the right plot show the harmonic approxi- 
mation to the trapping field. 



32 
3^ 



(39) 



The time evolution of the Gaussian wavepacket shows 
a strip of < 1^ ~ ro in which the density matrix 
decays much slower. This is a consequence of the C_ 
dependence of the decay rate shown in Fig. The spa- 
tial correlation length of the wavepacket can hence be 
read of Fig. ^ as Ic ~ tq. In addition, a damped os- 
cillation in the relative coordinate C_ is arising which is 
getting more pronounced for larger values of C,+ . Fig- 
ure jSl shows cuts along the C_ -direction for different val- 
ues of C+. The origin of the damped oscillations is the 
fc-dependence of the equilibration and decoherence mech- 
anism described by F^ in Eq. H34|) . To demonstrate the 



B. The double-wire trap 

The second configuration which we discuss is a double- 
wire trap 20] . A system of two parallel trapping wells is 
of special interest for high-precision interferometry pol 
I21I |2^ or beam-splitter geometries The effect of 

decoherence is one of the key issues in these experiments. 

Form and characteristic of the double- wire trap is de- 
scribed in [23, so that we will only briefly introduce 
the main features of the double-wire trapping potential. 
Figure ^ shows the setup for the trapping field, which 
is generated by two infinite wires running parallel to 
the z-axis, separated by a distance d, and a superposed 
homogeneous bias field -B^ias perpendicular to the cur- 
rent direction. A second bias field pointing along the z- 
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direction, B^'il^, is added in experimental setups to avoid spin flips. Thus the trapping field is 
I 



B(x) 



/io/ 
27r 



















^ bias 



(40) 



There are two different regimes for the positions of 
the trap-minima. Defining a critical wire separation 

yQ = ^qI / {2T:B^2,a) '-^^ distinguish the situa- 
tion of having d > 2yQ where the trap minima are lo- 
cated on a horizontal line at Xq^^' = ^\/(P /A — and 
yo = yoy a-iid the situation of d < 2yo where the trap 
minima are positioned on a vertical line at = and 
2/o^ = 2/0 i i/j/g — (P /4,. The two trap minima overlap 
and form a single minimum for d = 2yQ. 

Further on we will restrict ourselves to the regime d > 



2yQ since in this configuration we have two horizontally 
spaced, but otherwise identical trap minima. The two 
trap minima form a pair of parallel waveguides. This 
kind of geometry has been suggested for an interference 
device iH. Fi gure IHl shows a contour-plot of the 
double trap in the d > 2yo regime. 

We now analyze the double-wire setup along the lines 
of Section IIII Al Firstly, an explicit expression for the 
correlation function {5Sim{z)6Snj{z')) is needed. After 
some calculation we obtain 



{y - yo) - 



yaix-x'^) 



{Xq +£72) 

where the spatial correlation is now given by 



Xm(ri)Xn(r'_L) 



(y' - 2/0) - 



Vo(.x' - xg) 
[xq + e-yl) 



(41) 



dz 



^0-^72 ) +yl + ^^ 



2^0-^7 2 ) +yo + (2-2'-^) 



(42) 



Here, a is the trap label corresponding to i and to, and 
(3 the trap label corresponding to n and j, the indices 
which occur in {dSi„i{z)5Snj{z')). The wires are located 
at X — dj^/n — =F'^/2 and is defined as = — 1, 
= The result Eq. (|41|l uses the assumption that 
the extent of the wavefunction w is much smaller than y^ 
and much smaller than the trap minimum offset |a;Q^^| 

from the i-axis, i.e., w <C ja^o^^l w <^ y^. This 
approximation breaks down if the wire separation d ap- 
proaches the critical separation distance 2yo at which the 
two trap minima merge. 

Instead of going into details of the derivation we will 
discuss the different contributions qualitatively. As the 
system is invariant under mirror imaging at the y-z- 
plane, corresponding to an invariance under interchange 
of the left (L) trap and right (R) trap label, only the fol- 
lowing three contributions have to be distinguished (see 
Fig. 



(a) The influence of current fluctuations in the left wire 
onto an atom in a state which is localized in the left 
arm of the trapping potential (Fig.Et)- 

(b) The influence of current fluctuations in the left wire 
onto an atom which is in a superposition of a state 
localized in the left arm and a state localized in the 
right arm of the trapping potential (see Fig. Et")- 

(c) The influence of current fluctuations in the left wire 
onto an atom in a state which is localized in the 
right arm of the trapping potential (Fig. Et)- 

All other contributions are obtained by symmetry op- 
erations. Contribution (a) is given by the term propor- 
tional to J^^ and is the dominating decoherence source 
in the regime d > 2yo, as all other contributions {J^j^, 
Jrb. shown in the third and forth column of Table are 
suppressed by orders of yi^i/d. 
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FIG. 7: Contributions of the noise in the left wire to the 
decoherence of pij{z, z' ,t). The figures show only the coupling 
of one point on the wire to the atomic state. To get the total 
contribution to the decoherence an integration over all wire 
elements coupling to the atomic state has to be performed. 







JtuiC-) 






1 


2 

y„a 


1 


C- « yo 


1 


2 f -1 3Vo\ 


1 


yo < C- < 


4 


2 3CM 


1 


d< C- 


4 




4 



TABLE I: Estimate of the contributions to the decoherence. 
The columns for Jll, JlR) a-nd J^^ correspond to contribu- 
tions (a), (b) and (c) of Fig. Q respectively. The >/2fl(C-) a-re 
given by Eq. ^ as a function of — z — z' . The approxi- 
mation assumes yo <C d. 



The contribution in Fig. [Tja is of particular interest as 
it describes the cross-correlations between the noise in 
the left and right trapping well i.e. given by terms pro- 
portional to J^^ and Jf'fl. Current fluctuations in one 
of the wires give rise to magnetic field fluctuations in all 
the trapping wells. Magnetic field fluctuations are there- 
fore not uncorrelated. These correlations are however 
suppressed for yo <^ d as can be seen in Tabled 

Contributions of type (c), given by terms proportional 
to Jj^ji, are the smallest contributions to decoherence as 
they describe the influence of the left wire current noise 



onto the atomic cloud in the more distant right trap. 
Terms of type (c) are thus negligible for wire distances 
d ^ 2yQ as the magnetic field decreases with l/r and 
fluctuations in the trapping fleld are dominated by the 
nearest noise source i.e. contributions of type (a). In 
analogy to the proceeding Section Fill Al we are going to 
approximate the transversal wavefunctions Xn in Eq- l|4HI 
by 2d-H0 states. However one has to be careful as this 
approximation holds only if the two trapping wells are 
sufficiently separated such that the mutual distortion of 
magnetic field is small. If the distance between the two 
traps gets close to the critical value d — 2yo, a double- well 
structure arises and the harmonic approximation breaks 
down. We estimate the validity of the approximation by 
calculating the distance dc at which the local potential 
maximum, separating the two trapping wells, is of the 
order of the ground state energy of the isolated single 
trap. The trap frequency for the double- wire configura- 
tion extracted from the harmonic approximation is 



d d 



(43) 



where w is the single-wire trap frequency given in 
Eq. H26|) . Thus the ground state energy = Tilo of the 
single-wire trap corresponds to the limit of infinite sep- 
aration d — oo. This energy Eq is the maximum value 
for the ground state energy Tiujd of the double-wire trap 
under variation of d. Hence E^ is an upper bound for 
the ground state energy of the double well in the case 
d > yQ. Using the procedure described above, we get 
the condition d > VSyo, restricting the region in which 
the transversal potential can be approximated as two in- 
dependent 2d-HOs. The dashed line in Fig. shows an 
example for the harmonic approximation to the trapping 
potential of the double-wire setup. 

Performing the integration over the transverse coor- 
dinates rj^, using 2d-H0 states for Xn, gives transition 
matrix elements 



{dS.^iz)6SnAz')) = 4Ao^ ^ J2^(z - z') 



(44) 



(Xq -t 2) 



{Xq + e-y 2 ) 



where the functions Sij have been defined in Eq. (|28|l . 
For simplicity we choose the current in the left and right 
wire to be the same for the derivation of Eq. H45|l . 

Comparing the matrix elements obtained in Eq. H45|) 
with those in Eq. H27|) for the single-wire trap, we find 
the following differences. First of all it has to be noted 



that the transitions are no longer restricted to the y- 
components of the 2d-H0. This is a direct consequence 
of the geometry as the positions of the trap minima are 
now sensitive to current fluctuations in x- and y-direction 
whereas the x-position for the single-wire trap minimum 
was independent of the current strength in the conductor. 
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Secondly, the prefactor of the transition matrix element 
in Eq. (|45|l is now a function of the spatial coordinates 
Xq, taking the geometric changes of the trap minimum 
position under variation of the wire separation d into 
account. 




FIG. 8: Schematic setup of an interferometer for cold atoms 
using a double- wire trap. A wavepacket is coming in from the 
left single waveguide formed by a double wire trap with a wire 
spacing of d = 2yo. The wavepacket is split into a coherent 
superposition of states localized in the left and right arm of 
the double waveguide. Inducing a potential difference, e.g. 
by applying a current difference 51 in the left and right wire, 
gives rise to a phase shift. Remerging the wavepackets results 
into an interference pattern in the atom density. For a more 
detailed description of atom interferometers of this type see 
Ref. 113. 

Inserting Eq. H45|l in Eq. H21|) leads to a set of equations 
for the dynamics of the averaged density matrix. For 
simplicity we restrict ourselves to the ground state and 
first excited state of the 2d-H0. Having in mind a spatial 
interferometer as suggested in Ref. |2l|, the most inter- 
esting quantity is the density matrix element describing 
the coherent superposition of a state in the left arm with 
a state in the right arm. Figure |H| shows the geometry of 
the interferometer. A wavepacket is initially prepared on 
the left side in a single trap. The wavepacket propagates 
towards the right and is split into a coherent superposi- 
tion of a state in the left arm with a state in the right 
arm. Finally, applying a small offset between the cur- 
rents in the left and right wire gives rise to a phase shift 
between the two wavepackets which results, after recom- 
bination, in a fringe pattern of the longitudinal density. 
For a more detailed description of atom interferometers 
of this type we refer to ^21]. Loss of coherence between 
the wavepackets in the left and right arm, due to current 
fluctuations in the conductors, will decrease the visibility 
of the interference pattern. We will thus study the decay 
of the components of the density matrix off-diagonal in 
the trap index. We will suppress the trap label (L/R) 
in all following expressions, bearing in mind that the left 
(right) index of p will always denote a state in the left 
(right) trap. 

We will assume that the incoming wavepacket is ini- 
tially split into a symmetric superposition of the left and 
right arm. In addition we will choose the currents in the 
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FIG. 9: Decay time in the double- wire trap. T1/2 is defined 
as the time it takes until Pqq(z, z' ,t) has decayed to half of 
its initial value at t = 0. The graph shows its dependence 
on the wire separation d if all contributions in Eq. 1)46^ are 
taken into account (solid line), if the cross-correlations /3i are 
neglected (dotted line), and in the limit of d S> yo, Eq. (I52II 
(dashed line). The dash-dotted line is the ri/2-time using the 
decay rate Eq. (1381 for the single-wire configuration. Inset: 
r(t) = ln(poo)/t at d = 3yo- The trap-chip surface distance 
taken for the plot is yo = lO/im and all other trap parameters 
correspond to those used in Fig. 

left and right wire to be same, 1^ — Ir^ which does not 
restrict the validity of the obtained result for the deco- 
herence effects but keeps the equations simple as the full 
symmetry of the double- wire geometry is still conserved. 
The dynamics of p{z,z',t) is obtained from Eq. I|21|l . 
leading to a matrix equation in analogy to Eq. (|31|) but 
with a density matrix vector 

P-{p'o"o,pTi,pll,p'ilpl'i) , (45) 

where the upper index denotes the i-state and the lower 
index the y-state of the 2d-H0. The left (right) pair of 
indices refers to a state in the left (right) well. The matrix 
A is 



/ T,t=i 


Pi 


P2 


Pz 


Pi 


\ 




ai + as 





ag 






P2 





a2 + a4 


"5 


ag 






ae 




ai -|- a4 







V Pi 


"5 


as 





a2 + as 


/ 



(46) 

Here, terms labelled by ai are contributions of the type 
shown in Fig. plus the corresponding term of type 
Fig. 13:. Terms denoted by Pi are cross-correlations as 
described by contributions shown in Fig. [7|3. The terms 
abbreviated as a do not depend on the longitudinal coor- 
dinate z, z' whereas the cross-correlations /? are functions 
of {z-z'). 

We will assume that the extent of the wavepacket 
Wz along the trapping well is Wz d. As the cross- 
correlations Pi vary on a length scale of approximately 
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d, we can replace Pi{z — z') in Eq. (|46|l by the constant 

A(o) 0. 

Further on we will assume that a,t t — Q the density 
matrix is given by: 

p(i = 0) = (p!]SJ(t = 0),0,0,0,0) . (47) 

Figure El shows the time T1/2 it takes until p^{t = 0) has 
decayed to half of its initial amplitude as a function of 
the wire separation d. We assume the average currents in 
the wires to be constant and we vary only the separation 
length d. Thus the position of the trap minima given by 

2/0 = 2/0 , (48) 

(49) 



L/R. 



move horizontally to the surface, leading to a wire-to-trap 
distance of 



ro 



d 
71 



1-(1-4| 



(50) 



The solid line in Fig. |51 shows the T1/2 time obtained 
from the solution of Eq. H31() using Eq. H46|l in the approx- 
imation I3i{z — z') « /3i(0) for Wz ^ d. Cross-correlation 
terms (3i play however only a minor role for the 2\/2 time 
as can be seen from the dotted lines in Fig.|51 showing the 
case where the Pi terms are neglected in Eq. (|46|l . The 
cross-correlation terms, shown in Fig. \7}p, include terms 
describing positive correlation between potential fluctu- 
ations, suppressing the decoherence, but also negative 
correlation between potential fluctuations which enhance 
the decoherence. Figure ^| shows schematically the po- 
tential fluctuations 6V induced by a current fluctuation 
SI in the left wire. 

The inset of Fig. |51 shows the change of the decay 
rates over time. Taking the full matrix Eq. (|46|l into ac- 
count (solid line) there is a change of decay rate Ta+pit) 
towards a slower decay rate for t '3> r~j|^^(0). If the 
cross-correlations are ignored, the decay rate is constant 
Tq « ra+p{0)- Inspection of A, given by Eq. (|^ . shows 
that neglecting the terms f3i decouples from the ex- 
cited states leading to an exponential decay with a single 
decay rate. 

We finally discuss the decoherence for large distances 
d between the wires, i.e., d ^ yo- For d ^ yo it is suf- 
ficient to take only the influence of 6Il {SIr) onto the 
left (right) trap (given by contributions of type Fig. E^) 
into account, as all other contributions are strongly sup- 
pressed by orders of t/o/d (see Tabled). Hence the matrix 
Eq. 1)46(1 decouples Pqq from the excited states and the re- 
maining equation for Pqq can be solved analytically. The 
solution is 



-rt. 



rHt MO 



p°°(z,z',0), (51) 



where H is the free Hamiltonian as defined in Section 
III Al The decay rate is given as 



Stt \xo\' 



1 



{\xo\+d/2Y 



1 



(52) 







1+ 51 

FIG. 10: Influence of current fluctuations 51 in the left wire 
onto the potential of the trapping wells V{x, y) = (S'|/j.|5')B. 
The potential fluctuations 5V{x,yo) — 5V{—x,yo) are reflec- 
tion symmetric about the j/-z-plane for \x — Xq^^\ ^ yo- 
Hence the fluctuations in £-direction are conserving the sym- 
metry of the trapping potential. This leads to positive correla- 
tions in {5Sim[z)5Snj(z')) suppressing the decoherence. The 
fluctuation in y-direction, 5V{x^ ,y), SV{x§,y), show nega- 
tive correlations giving rise to an increase in the decoherence 
rate. 



In the limit d^ yo, Eq. (|52|l reduces to the decay rate ob- 
tained for the single- wire configuration, Eq.ljJSJ, shown 
by the dash-dotted horizontal line in Fig. |51 However, 
there is no spatial dependence of F along the longitu- 
dinal direction (z — z') in the case of widely separated 
trapping wells. Plotting the T1/2 time for the approxi- 
mation given by Eq. (|51|1 (dashed line in Fig. O using 
the decay rate Eq. H52|) . shows a good agreement of the 
approximation with the exact solution as soon as d is of 
the order of several yo. Comparing all three graphs of 
Fig. 1^ one can conclude, that the increase in with 
decreasing separation d is due to geometric changes of 
the trap positions. An increase of the wire-to-trap dis- 
tance ro due to a decrease of d, leads to an increase in 
the Ti/2 time as the rate of the dominating decoherence 
source scales with 1/rg. 



IV. CONCLUSION 

Using the density-matrix formalism, we derived an 
equation of motion Eq. ifT^ to describe the consequences 
of current fluctuations on a cold atomic cloud in a micro- 
chip trap. The model allows the description of deco- 
herence in multiple wire traps |3lL Is^ , as well as more 
complex guiding systems required for the implementa- 
tion of a beam splitter [l^ or interference experiments 

[loliillil. 

The atom trap was modelled as a multi-Channel Id- 
waveguide, where different channels describe different 
transverse modes of the waveguide. Assuming Id-wires 
on the atomic chip, we examined the influence of cur- 
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rent fluctuations along the wire onto the coherence of the 
atom cloud in the waveguide. We found that important 
contributions to the decoherence arise from transitions 
to neighbouring excited states. 

We used this model to examine decoherence effects 
for two specific trapping configurations: the single-wire 
waveguide and the double-wire waveguide. In both con- 
figurations, decoherence of the ground state was dis- 
cussed, taking processes from transitions to the first ex- 
cited states into account. 

The single-wire trap showed for the ground state a de- 
coherence rate T which scales with the wire-to-trap dis- 
tance ro as I/tq. The potential fluctuations are corre- 
lated over a length scale rg. As a consequence the de- 
coherence rate F is a function of the relative coordinate 
z — z' . Using trap parameters based on present experi- 
ments js^ l we obtained decoherence rates for ro = 10/ini 
of the order of T w 0.03s-^ 

Extending the system to a double-wire waveguide en- 
abled us to study decoherence for atoms in a superpo- 
sition of a state localized in the left arm and a state 
localized in the right arm. This superposition is the ba- 
sic ingredient for interference experiments. Approaching 
the two trapping wells, as is necessary for the splitting 
and merging of the wavepackets, showed a decrease in 
the decoherence rate. The decrease arises mostly due to 
geometrical rearrangements of the trap minima in the 
system. Cross-correlation effects proved to be of minor 
importance in our model. We found an explicit expres- 
sion for the decoherence rate Eq. (|52|l in the limit of a 
wire-to-trap distance ro much smaller than the separation 
of the two wires d. This decoherence rate approaches the 
value for the single wire waveguide, T = A{0), for de- 
creasing ro/d. 

The decay rates extracted in this model are small for 
distances realized in present experiments [3^ . However 
we expect that further improvement in micro-structure 
fabrication and trapping techniques will decrease the 
trap-surface distance to scales at which the decoherence 
from transitions to transverse states may have a consid- 
erable influence on trapped atomic clouds. 
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APPENDIX A: DERIVATION OF THE 
PROJECTED POTENTIAL CORRELATOR 

This section gives a derivation of the correlation func- 
tion {5Sim{z)5Snj{z')) ■ A general expression of the po- 
tential fluctuation correlator (5F(x)(5t^(x')) for N par- 
allel traps will be derived, from which {5Sim{z)5Snj{z')) 
can be calculated using Eq. Finally Eqs. - t^ 

and Eq. H41I) , Eq. 1)42(1 can be obtained by specifying the 
geometries specializing to the single and double- wire trap 
configurations. 

The fluctuations of the trapping potential 5V is in- 
duced by current noise in the wires which gives rise to a 
fluctuating field (5B. The trapping potential 5V and the 
magnetic trapping field (5B are linked by Eq. which 
enables us to express {5V{:x)5V{'x')) in terms of magnetic 
field fluctuations: 

(<5l^(x)5F(x'))-^(5(x)|^,|5(x)) 

X {S{^)\^^,\S{^)){5B,{^)5B,{^')) . 

(Al) 

Thus the first step is to find the expression for 
{6Bi{x.)SBj{:x.')) as a function of the current noise. We 
start by calculating SA and using 

M.(x) =M/d3igd|, (A2) 
<5B.(x) ^J2Jk^^Jk£-SAk{^). (A3) 

The current density in the set of Id- wires is 

5j(x) = J2 SI^i^)S{x - d^)S{y)z , (A4) 
7 

where d^ denotes the i-position of the 7-th wire. The 
sum over 7 reflects the fact that the total current density j 
is a sum of M contributions to the current density arising 
from the set of M wires. We will now insert Eq. I|A4|I in 
Eq. (|A2p and Eq. l|A3p to obtain an expression for the 
magnetic field induced by the current in the wire array: 

^ V / (A5) 

J [{x - d^f + + {z - zff'^ 

The correlation function for the current density is now 
assumed to be of the form 

{5Io,{z)5Ifi{z')) = AkYiT,s(rA5{z - z')5o,p (A6) 

for the reasons discussed before in Section III CI Using 
Eq. (|A6|) in combination with equations I|A1|I and (|A5|) 

gives the desired relation for the correlation function of 
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the potential fluctuations 

{SV{^)SV{^')) = 4fcBTeffaA 

X ^^(5(x)|/.,|5(x))(5(x')|a*,|5(x')> (A7) 

7 »i 

X yJ(x,x')J^(x,x'), 
where the following abbreviations have been introduced 



/ 



y^(x,x') = 



J^(x,x') 



yy 



-{x — dj)y' (x — d^)(x' ■ 




d^) 



(A8) 



dz \{x — d^Y + y"^ + ~ z) 



[{x' -d^f + y'^ + {z' -if 



(A9) 



Equation (|A9|) can be further simplified if the transver- 
sal positions x^'^ and y''^ are replaced by the position 
of the trap minimum x'^ and where a is the trap la- 
bel denoting the trap in which the wavefunction is local- 
ized. The replacement of the transversal coordinates by 
its trap minima positions is a good approximation as the 
transversal widths w of the trapped atomic clouds are in 
general much smaller than the wire-to-trap distance tq. 
We have thus reduced J'' to a function which now only 
depends on the difference C_ — z — z': 



.,2 



J2«(C-)= / d~z' {x^,-d,f + y^'+~z' 



(AlO) 



X [xl^Q - d. ^2 



)^+yr + (c--gr 



where we shifted the integration variable to z' = z — z. 
Applying the formula to the single and double wire con- 
figurations leads to Eq. H25|l and Eq. H42|l respectively. 

To obtain {5V(-x)5V{-x')) using Eq. I^Kl\. we still need 
to calculate the mean value of the atomic magnetic mo- 
ment (S'(x)|/x|S'(x)). We assume that the magnetic mo- 
ment follows the magnetic trapping field adiabatically 
which is reasonable as long as the Larmor precession 
— I^-bB/H is fast compared to the trap frequency uj. 
Calculating the spinor |>S'(x)) for an atom having spin 
F — 2, either by considering the small corrections of 
the transversal magnetic field to i^^^as perturbatively or 
by calculating the rotation of the spinor as the atomic 
moment follows the trapping field adiabatically, results 
in the following expression for the spatial dependence of 
|S'(x)) for small deviations from the trap minimum: 



|^(x)> = |2,2>-H^^iM±lM^|2,l) 



r(^) 



(All) 



The spin states are denoted as \F,mF) and the spin quan- 
tisation axis is chosen along the z-axis. Making use of 
Eq. HA11|) leads to 



(5(x)|/x|5(x)) = 2/iBgF 



B(x) 



(^) 



(A12) 



which is an approximation to the mean value of the 
magnetic moment to first order in B±/B^^^}^^. Inserting 
the expression for the magnetic moment Eq. ljA12|) into 
Eq. ljA7|l and inserting the specific magnetic trapping 
fields for the single- wire trap, Eq. H22|l . or the double- 
wire trap, Eq. H4()(l . we finally obtain the results for 
{6Sim{z)SSnj{z')) given by Eq. for the single- wire 
setup and by Eq. H41|l for the double-wire setup. 
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